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Summary 

The stability/instability condition of a turbine rotor with axisymmetric supports is determined in the 
presence of gyroscopic loads and rub-induced destabilizing forces. A modal representation of the turbine 
engine is used, with one mode in each of the vertical and horizontal planes. The use of nonspinning rotor 
modes permits an explicit treatment of gyroscopic effects. The two linearized modal equations of motion 
of a rotor with axisymmetric supports are reduced to a single equation in a complex variable. The result- 
ing eigenvalues yield explicit expressions at the stability boundary, for the whirl frequency as well as the 
required damping for stability in the presence of the available rub-induced destabilization. Conversely, 
the allowable destabilization in the presence of the available damping is also given. 


Introduction 

Simulation of blade to case rubbing is critical to the complete understanding of the dynamic behavior 
of rotordynamic systems such as aircraft engines. Aircraft engines experiencing unbalance and maneuver 
loadings, blade-outs, or mechanical misalignments may incur blade to case rubbing which may adversely 
affect the dynamic behavior of the engine and supporting structures. In some situations, the rubbing may 
become excessive and may even lead to structural instability and in severe cases catastrophic failure of 
the structure. 

The purpose of the present paper is to determine the conditions for stability/instability of a turbine 
engine with rub-induced forces in the presence of gyroscopic forces. To develop a direct determination of 
the stability/instability conditions, the equations of motion with rub forces obtained in reference 1 are 
used as the starting point for the present investigation. In this study, the vertical and horizontal planes are 
each represented with one nonspinning vibration mode. Stability/instability is then investigated by exam- 
ining the corresponding eigenvalues of the characteristic equation. 

The thrust of this work is with a rotor with axisymmetric support stiffness. This model allows for a 
complex- variable representation of the rotor system, and a more direct determination of conditions at the 
stability-instability boundary. The notion of neutral stability at the boundary permits obtaining two real 
equations from the complex eigenvalue problem. Equations are derived for the possible whirl frequencies 
and the required damping for stability (or the allowable de-stabilization), at the stability/instability 
boundary. 
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Blade Rub Forces and Moments 


Limited data and experience with rub tests on turbofan engines and test rigs indicate that blade tip rub 
forces, resulting from blade tip and case interactions, are in a direction that is perpendicular to the blade 
tip chordline. This phenomenon is similar to the behavior of an airfoil traveling through a fluid where the 
airfoil motion generates a lift and drag forces. From tests, the tangential rub blade force was found to have 
an axial component that would move the blade forward (or aft), and a circumferential component that 
opposes the rotor rotation. The hydrodynamic analogy expresses the blade rub force as proportional to the 
blade tip’s mean velocity and the instantaneous angle of attack or incidence. It differs from the typically 
used friction model of rub in that the friction model only generates a force in the direction opposing the 
rotor spin (i.e., drag) while the plowing model generates this drag force, as well as a “lifting” force that is 
perpendicular to the direction of rotor spin. A more detailed description of this rub model is provided in 
reference 1. 

An analytical expression, analogous to the dynamic pressure force from a change in momentum is 
used to develop the plowing model. This model is consistent with expressions of impact forces and 
hydrodynamic pressure. This was necessitated by the observed directionality of the blade tip “tangential” 
rub forces (i.e., perpendicular to the blade chordline). 

The hydrodynamic force is defined as: 


f = Dy 2 R 2 y 0 A=DV 2 y 0 A (1) 

where / is the total tangential force on the blade from rubbing, vj/ is the rotor rotational speed, R is the 
blade tip radius, V is the blade tip velocity, y 0 is the blade chordline angle of incidence in the rotating 

frame, and A is the frontal area of the rubbing material (i.e., chord length x rub depth). D is a coefficient 
that must be determined from test or some other analysis procedure. It is most likely a function of several 
factors such as the material’s Young and Shear moduli, density, hardness, machinability, feed depth, 
dynamic shear strength, smoothness, and probably others. Figure 1 provides a schematic of the rotor blade 
and related variables. 

The loads in the fixed coordinates system are: 

f x = f sin y o sin (I - K r x (2a) 

f y = ~f sin y (, cos P - K r y (2b) 

fa = -f cos Yo cospfi (3a) 

./ 0 = f cosy q sinfiR (3b) 

Where Kr is the radial stiffness. Noting that: 

sin B = v and cos B = -y 
r /r r /r 

Where r is the blade radial displacement. Substituting into equations (2) and (3), the rub loads 
become, 


f x = Dy 2 R 2 y 0 ysiny 0 - K r x 


(4a) 
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(4b) 


f y =-D\\i 2 R 2 y 0 xsiny 0 -K r y 

fa = -D<\r 2 R 3 y 0 xcosy 0 (5a) 

fQ = D\\i 2 R 3 y 0 ycosy 0 (5b) 


Modal Transformation of the Rub Loads 

As a first step to simplifying the equations of motion, the translational and rotational rub forces in 
equations (4) and (5) are transformed to modal coordinates. The following considerations are used to 
obtain the modal rub forces: 


1. There is one modal function in each of the horizontal ( x , a) and vertical (y, 9) planes. 

2. The rotational (a, 9) deflections are the slopes of the bending displacements. 

3. Because of the (x, y, a, 9) sign conventions, it is seen that: 


9 = 


dy 

dz 


and 


dx 

a = + — 

dz 


4. The modal functions in the fixed axes are nonrotating. 

5. The modal functions are mass normalized. 


Now let the modal displacements be as follows: 


x = xh{t) 

(6a) 

y = yv(f) 

(6b) 

a = ah{t ) 

(7a) 

CD 

II 

1 

CDI 

< 

(7b) 


Where x , y , a , 9 are mode shapes and h (t) and v(t) are horizontal and vertical plane displacement 
generalized coordinates, respectively. 

Following the Principle of Virtual Work, the modal forces in the horizontal and vertical planes are: 


fh=j{ x fx + a fa) ds 

(8a) 

fv = Jfe -We)ds 

(8b) 


Where the rub forces are integrated over the tangential extent of the rub, “s” then summed over all 
rubbing blades. 
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Making the physical-to-modal displacement substitutions, the modal rub forces become: 


fh= j{0¥ 2 fl 2 Yo(*y sin Y(p ; (t)- Rolx cosyQh(t))~ K,-x 2 h(t)jds 

f v = J D\y 2 R 2 yo{xy sinjoh{t)+ RQycosyQv(t^j-K r y 2 v(t)^ds 
Defining the following modal integrals: 

K x - j K r x 2 ds 
S .xy - J xyds 
s$y — J Qyds 

— f 2 

Ky - I K r y ds 

S ax = J axds 

Then substituting into equations (9a) and (9b) yields: 

fh - D <r Rl yo Krv sin y 0 v(f)- RSax cosy 0 h(t)}~ K x h(t) 

fv = ■ d ¥ 2 ^ 2 Yo {“ s xy siny 0 h(t)-RS Qy co,vy 0 v(f)}- K y v(t) 

The following variables are defined to ease the required algebraic manipulations: 

5 0 = D\\i 2 R 2 y ()S xy sin y () 

Sh = D\y 2 R 3 y 0 S ay cos y 0 
S v =D\y 2 R 3 Y 0 S Qy cosy 0 

Substituting into equations (10a) and (10b) the modal rub force equations are reduced to: 

fh =S 0 v(t)-K x h(t)-S h h(t) 
fv =-S 0 h(t)-K y v(t)-S v v(t) 


(9a) 

(9b) 


(10a) 

(10b) 


(11a) 

(lib) 
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The Modal Equations of Motion 


The modal equations of motion are obtained by combining the modal rub forces (eqs. (11a) and 
(lib)) with the modal inertia, stiffness, damping, and gyroscopic forces. The latter forces are obtained by 
using the translational and rotational values of the horizontal and vertical modal vectors to transform the 
physical inertia, stiffness, damping, and gyroscopic properties to modal coordinates. Performing these 
transformations and using the modal rub forces derived in the previous section yields the following modal 
equations of motion (see also refs. 2 and 3): 

h + c ojh + 2co ,& h h + 2\j igv - f h (12a) 

v + a>vV + 2co v ^ v i> - 2ygh = f v (12b) 

Substituting the rub forces and collecting coefficients of the same dependent variables we have the 
following: 


h + (to l + K x + S h )h + 2co h % h h + 2vj igv - S Q v = 0 (13a) 

v + (co; + K y + S v )v + 2co v ^ v v - 2\\igh + S Q h = 0 (13b) 

In general, a large portion of rotors have nonaxisymmetric supports so the modal stiffness and damp- 
ing differs in the horizontal and vertical directions. If the horizontal and vertical modal displacements are 
not equal then the direct rub force coefficients S/„ and S v , K x and K y , are also not equal. However, many 

incidents of turbine engine whirling exhibit circular or almost circular orbits, even for nonaxisymmetric 
rotors. Furthermore, one may also be able to construct highly elliptical orbits from identical uncoupled 
horizontal and vertical mode shapes even when the corresponding modal stiffnesses or frequencies are 
grossly different. Therefore, to simplify the modal equations of motion it will be assumed that the hori- 
zontal and vertical mode shapes are the same. However, the support stiffnesses and modal frequencies 
may be unequal. 

With this consideration the following additional simplification is made: 

S h = S v = Si 
K x =K y =K 

Where the local case and rotor radial stiffness is assumed to be isotropic and axisymmetric. Addition- 
ally, the modal damping is assumed to be identical in both planes. The damping symmetry is justified 
since damping values are difficult to obtain and in most analyses is an educated guess, at best. Assuming 
symmetric damping: 


Cj 2oi/,qy, 2co v ^ v 


Finally, the modal equations of motion are: 


h + (eo/t + K + >S) j /7 + t^h + 2\j igv — 5qv — 0 


(14a) 
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v + + K + S\ j v + Cy - 2\j igh + SqH = 0 


(14b) 


Where the support nonaxisymmetry is reflected in the horizontal and vertical modal frequencies, ©/, 
and © v . 


Stability of Rotor With Axisymmetric Supports 

It had been shown by various investigators such as Gallardo (ref. 4), Smith (ref. 5), Ehrich (ref. 6), 
Alford (ref. 7), et al., that the least stable rotating system is axisymmetric. Thus one may initially perform 
a stability analysis of a rotor with an axisymmetric support. Then, if the axisymmetric system is stable, no 
additional analysis may be required, since a nonaxisymmetric supported system will always be more sta- 
ble than the identical rotor with symmetric supports. 

In the special problem of an axisymmetric rotor system the stability/instability conditions are obtained 
with greater simplicity and directness than for a rotor with nonaxisymmetric supports. The modal equa- 
tions of motion of a rotor with axisymmetric support stiffness are obtained from the previous section by 
making the modal frequencies equal. Thus: 


h + |coq + jSl + Sj j/z + Cfi + 2\j fgv — Sqv — 0 

(15a) 

v + |coq + jST + ^jv + Cy — 2\j/g/i + SqH — 0 

(15b) 


Where © 0 = ©/ 7 = © v . 

Equations (15a) and (15b) are combined into a single equation by introducing the complex variable z, 
where z = h + iv. Multiplying the second equation of motion (eq. (15b)) by i and adding it to the first 
equation of motion (eq. (15a)), and using the complex variable definition, we have 


z + (©o + K + Si )z + C,z - i2\\igz + iS 0 z = 0 
Assuming a solution in the complex exponential form: 

i'Kl i(A+iB)t -Bt iAt 
Z = Zoe =ZQe y ’ = zq? e 


(16) 


Where the coefficient z 0 and the eigenvalue X are complex quantities and A and B are real numbers. 
Furthermore: 


Real (X) = A = frequency (real quantity) 

Imag (X) = B = attenuation or amplification factor (real quantity) and: 
B > 0, stable 
B < 0, unstable 
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Stability/Instability Boundary of the Axisymmetric System 

Substituting the exponential solution into the complex variable equation of motion (eq. (16)), one 
obtains: 

|^ 2 -2^vi/g-(^+Z + 5 1 )^-/(^ + 5 0 )| z = 0 (17) 

The condition for the existence of a nontrivial solution requires that the polynomial coefficient of z 
must vanish: 

F(x) = j^A 2 - 2X\j/g -(coo + K + sJj- i(X£ + S 0 ) j = 0 (18) 

At the stability-instability boundary, the solution is simple harmonic, which implies that X = Real at 
the boundary. Therefore, since all the parameters in the Eigen polynomial (eq. (18)), are real and nonzero, 
then the following holds at the stability/instability boundary: 

Real part of F (X) = 0 
Imaginary part of F (X) = 0 

And we have two equations at the stability boundary: 

A, 2 -2A,\jfg-(<D§+* + 5^ = 0 (19a) 

XL, + Sq = 0 (19b) 

Equation (19a) is solved for X which gives the possible frequencies at the boundary. 


^ - W8 - V V 2 ^ 2 + (®o + K + 5i ) 

Using the frequency, X determined above and equation (19b), the damping required for stability is: 

X 

Since the term So is always destabilizing (S 0 > 0) and for a system to be stable, the system damping C 
must be positive (£ > 0). Therefore, X = X.< 0 at the boundary, because B = 0 at the boundary, and C, and 
So are positive. Furthermore, X - X+ > 0 is always stable because if/. > 0 then B> 0. 

From the above, the results for the rotor-rub stability/instability boundary are: 



Whirl frequency: X_ = vj/g - 



+ (coq + K + 



And the required damping, C, RE q 


So 
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Comparison of the available damping to the required damping (in the presence of a destabilization) 
gives a measure of the stability of the system: 

C AVAILABLE > C REQ > Stable 


C AVAILABLE ( ^REQ > unstable 


Recall that 5 0 and 5, contain the plowing rub coefficients, comprised of the rub constant, rotor speed, 
blade tip angle, and tip radius and Y , is the radial-rub restoring or snubbing force, which may or may not 
be dependent on the rub parameters. Recognizing: 


cot Yo = 


cos Yo 
sin yo 


Si may be written in terms of So as: 


S] = S Q Rp, where p ■■ 


cot Yo 


J xy 


Making the above substitutions into the expressions for the whirl frequency and required damping, we 
have, at the stability instability boundary, the following results. 


Whirl frequency: k_ = vj/g - -y vj i 1 g 2 + ^coq + K + S 0 Rpj (20a) 

£ 

Required damping: ^ RE q = 0 (20b) 

VS - “V vV + (coo + K + S 0 Rp ) 

The rotor shown in figure 2 is used to demonstrate the use of equation (20) for determining the 
whirl frequency and required damping for a known quantity of rub. The rotor consists of a massless 
cantilever shaft with a disk attached at the free end. The properties of the rotor are shown in the figure. 

The nonrotating frequencies and modes shapes are computed from a finite element analysis. The first 
modal frequency and mode shape is: 


9 

co 0 = CD/, = co v =3.071x 10 


rad / sec 


\a/ \e/ \ 7. 144 x 10 -2 / 

Substituting the above values, and the values for blade angle, disk radius, rub stiffness, and length, 
into equations (20a) and (20b), the whirl frequency and required damping are: 

X =0. 2644 y±ijo. 0699 \j/ 2 + 14.45 xlO 4 +37.04 Dvj/ 2 

_ 109.6£>vj/ 2 
^REQ~ r 
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Section A-A 

CD-04-82602 


Figure 1 . — Schematic of rotor blade, degrees of freedom and tip penetration. 


Disk 

M = 100 Ibm 
Ip = 20,000 lbm-in. 2 
Gravity = 386 in./sec 2 
Radius = 5 in. 

Blade angle = 30 deg 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

◄4 

Rotation 
vector 



Rotor shaft 
E = 30x1 0 6 lb/in. 2 
I = 20 in. 4 
- L = 40 in. 



Rub spring 
K = 3x1 0 3 lb/in. 
Rub length = 5 in. 


Figure 2. — Cantilever rotor. 


Figure 3 shows the result of plotting the whirl frequency and required damping as a function of rotor 
speed for a variety of dynamic rub coefficients. The whirl frequencies associated with both the forward 
(X+) and the backward (X_) mode are shown in figure 3(a), however, only the backward whirl frequency 
is significant since the backward whirl frequency is used to calculate the required damping for a stable 


NASA/TM— 2004-2 12974 


9 






system. The required damping as a function of rotor speed for various levels of rub is shown in 
figure 3(b). As expected and depicted in the figure, the required damping increases with both rotor 
speed and rub force. 


Allowable Destabilization for an Available Damping 

The damping relation may be solved for the allowable destabilization So, in terms of the available 
damping. In this context, equation (20b) is solved (with the available damping assumed to be known), for 
the destabilization that can be allowed that will make the initially stable system approach neutral stability 
or the stability /instability boundary. 

Thus the corresponding relation is: 


S 0 = -K 


And the converse expression for the allowable destabilization in the presence of the available damp- 
ing becomes, 


5 0 Movable Available jvS “ 'V + (®0 + * + S 0 Allowable R p) j ( 2 1 ) 

Equation (21) is complicated since S 0 is on both sides of the equal (=) sign. However, with some 
algebraic manipulation, the allowable destabilization can be obtained explicitly. Rearranging equa- 
tion (21) yields: 


0 Allowable 




( 22 ) 


Noting that So occurs skew-symmetrically in the modal equations of motion (eq. (13)), and that it has 
a positive sign (+) in the vertical modal equation and a negative (-) sign in the horizontal. Thus, in the 
mathematical context, only a positive (+) for So is admissible. The signs of So in the equations of motion 
reflect the unstable whirl direction, such that a positive sign in the vertical DOF (negative in the horizon- 
tal DOF) may indicate a backward whirl, while opposite signs would denote a forward whirl. This effect 
is also seen in the expression for the required damping. 

The assessment of system stability or instability can then be made by comparing the available destabi- 
lization to its allowable value, in the presence of the available damping. Thus: 

Sr. (So , stable 

w Available w Allowable 


So 


Available 


)S 0 


Allowable 


, unstable 
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Required damping Whirl frequency, Hz 




0 50 100 150 


Rotor speed, Hz 

Figure 3. — Rotor stability, (a) Whirl frequency at stability/instability boundary, 
(b) Required damping at stability/instability boundary. 
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Concluding Remarks 


The linearized modal equations of motion of a rotor with nonaxisymmetric stiffnesses with both gyro- 
scopic and rub forces have been obtained. The equations of the special case of a rotor with axisymmetric 
support have also been derived. The latter equations are used to determine the influence of the rub forces 
on rotor stability/instability boundaries. 

The two-mode representation (one in each vertical and horizontal plane) of the equations of motion 
was considered. For the axisymmetric support case, the modal equations of motion were reduced to one 
equation using a complex variable to describe the rotor motions. A method was derived to investigate the 
conditions at the stability/instability boundary. 

At the boundary, the whirl frequencies and required damping for the available value of the destabiliz- 
ing coefficient are obtained. The converse, which is the allowable destabilizing coefficients in terms of 
the available damping, is also derived. These expressions are in closed form, and provide a convenient 
way to assess the stability/instability implications of both rub force parameters and the rotor system’s 
modal inertial, damping, elastic and gyroscopic properties. 

The analysis for the axisymmetric rotor system may be considered as a negative stability/instability 
criterion, in that an axisymmetric system is the least stable configuration. Thus, if this analysis indicates a 
stable rotor, an analysis of a rotor with nonaxisymmetric support stiffnesses may not be required, since 
the latter configuration is always more stable than the former. 
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